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Abstract. A strategy for the systematic design of polymeric superlattices with tailor-made 
mini-bandgaps and carrier mini-effective masses is described and computationally implemented 
by means of an envelope crystalline-orbital method, which is a straightforward adaptation for 
molecules of the envelope-function approximation widely used in solid-state physics. Such 
strategy relies on the construction of 7r-conjugated periodic block copolymers from well- 
characterized parent polymers, in such a way that the above-mentioned electronic parameters 
can be predicted from the lengths of the blocks. Illustrative calculations for prototypical (PPP^- 
PDA y ) n superlattices demonstrate the plausibility of the strategy and the advantages of the 
computational implementation employed. 

To appear in International Journal of Quantum Chemistry 
1. Introduction 

A solid-state superlattice is a periodic arrangement of many nanometric layers of two, or more, 
different semiconductor materials [TJ [2J . The ionization energies and electron affinities of these 
materials are such that the effective potential-energy profiles for electrons and holes along the 
direction of growth take the forms of periodic arrays of quantum wells. The barriers between the 
wells are thin enough so that their otherwise isolated single-particle levels become resonantly 
coupled by quantum tunneling. Such coupling mixes these levels, giving rise to so-called energy 
mini-bands, in an analogous fashion to the appearance of energy bands in a crystal. In this 
sense, a superlattice behaves as an artificial crystal with a much larger period (of the order of 
10 nm) than any real crystal, and, consequently, with a much smaller (mini-)Brillouin zone in 
reciprocal space. The mini-bandwidths and, consequently, the carrier mini-effective masses are 
determined by the thickness of the barriers, whereas the mini-bandgaps are determined by the 
depths and widths of the wells. Therefore, the mini-band structure of a superlattice can be 
engineered by controlling these parameters [3]. 

Crystalline heterostructures such as this can be manufactured from inorganic [1] and organic [2] 
materials, employing sophisticated, and costly, epitaxial crystal-growth methods, like molecular- 
beam epitaxy and chemical or physical vapor deposition techniques [3]. These methods afford 
so good a control over the fabrication that the miniband structures can be tailored for specific 
electronic and optoelectronic applications [I]. 



This report deals with a different kind of superlattice, which is constituted by a quasi-linear 
7r-conjugated alternating block copolymer [5]. Now, the effective potential-energy profiles for 
carriers assume the forms of periodic arrays of tunnel-coupled quantum dots. (To conform to 
the current usage in the semiconductor literature, it is more appropriate to employ the term 
'quantum dot' than 'quantum well' in this case, since in the former the carriers are confined in 
the three directions, whereas in the latter they are confined in one direction only.) Hence, this 
organic heterostructure can be considered as a molecular onedimensional analog of the afore- 
mentioned solid-state layered superlattice. 

Molecular organic heterostructures such as this present the following advantages over their solid- 
state counterparts: First, in principle, it is possible to construct a virtually unlimited variety 
of them, due, in turn, to the virtually unlimited number of monomers available, whereas the 
number of parent crystalline materials available is limited [3]. Second, the "interface" between 
two blocks is a chemical carbon-carbon bond, in contrast to the interface between two crys- 
talline layers where it is required that the materials exhibit similar lattice constants [3], which 
notably limits the variety of heterostructures that can be obtained. Third, organic copolymers 
[6] and heterostructures [7] can be obtained through chemical-synthesis methodologies, whose 
implementation is easier and cheaper than the combination of epitaxial crystal-growth and litho- 
graphic methods commonly used for the manufacture of crystalline heterostructures. 

In the light of these observations, the general objective of this contribution is to point out 
that molecular quantum-dot superlattices, constituted by quasi-linear 7r-conjugated alternating 
block copolymers, afford an attractive alternative to solid-state quantum-well superlattices for 
the design of materials with engineered bandgaps and conductivities. Specifically, it shall be 
shown that, by appropriately selecting the parent polymers and defining the block lengths, the 
mini-bandgaps and carrier mini-effective masses of these molecular organic heterostructures can 
be tailored at will. Hence, this strategy is complementary to the one presented in the previous 
contribution to this issue [8], where the blocks are so short that the copolymer does not behave 
like a superlattice but instead like a new polymer, whose bandgap and effective masses can 
be controlled by the mole fractions of the parent monomers. (It should be clarified that these 
regular copolymers have also been called 'superlattices' by other authors [9]; however, this is 
not recommendable since these systems do not exhibit mini-bands). 

The implementation of this design strategy will be purely computational. In particular, firstly 
the parent polymers will be characterized by means of quantum-chemical electronic structure and 
bandstructure methods. Secondly, once the block copolymers are designed and assembled, their 
mini-bandstructures will be determined employing an envelope crystalline-orbital approximation 
ITTI [T2] , which is conceptually more attractive and computationally more economical for this 
task than any of the conventional atomistic quantum-chemical methods employed thus far for 
copolymers [SJE]. This method is a straightforward adaptation for molecules of the envelope- 
function, or effective-mass, approximation widely employed in solid-state physics [3, 13J. 

2. Design Strategy 

The design strategy consists of the following steps: (1) Choice of a pair A-B of parent polymers 
with an appropriate alignment of their valence-band tops, or highest unoccupied crystal orbitals 
(HOCO's), and conduction-band bottoms, or lowest unoccupied crystal orbitals (LUCO's). (2) 
Design of the repeat unit (quantum dot) of the superlattice, to guarantee that it possess the 
desired number of bound levels and range of separations between them. (3) Assembly of the 
superlattice, or alternating block copolymer, from the repeat unit of the previous step. (4) 
Determination of the mini-bandgaps and electron and hole mini-effective masses of the super- 



lattice, as functions of its structural parameters. (5) Construction of calibration curves of such 
quantities as functions of these parameters, from which the same quantities for an arbitrary 
superlattice can be extracted. The meanings of these steps are now explained. 

To illustrate the strategy, for step (1) poly-(p-phenylene) (PPP) and poly-diacetylene (PDA) 
were selected as parent polymers for the construction of quasi-linear block copolymers. Fig. [1] 
displays the HOCO and LUCO energies of these polymers, calculated using the methodology 
indicated in Section^ This alignment of the frontier CO's indicates that 'type F heterostructures 
[31 [J3] can be designed from these parent polymers (see below) . 
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Figure 1. Alignment of the PDA and PPP frontier CO's (band edges). Solid lines: HOCO's (valence-band 
tops). Dashed lines: LUCO's (conduction-band bottoms). 



For step (2), the three-block molecular heterostructure PPP/PDA/PPP is chosen as repeat unit 
for the superlattice. According to the envelope molecular-orbital (EMO) theory [51[TU l [T H [r2 | [r5]. 
the heteroj unctions between blocks can be considered as abrupt, so that the effective potential- 
energy profiles along the longitudinal direction for electrons and holes within the molecular 
frame take the forms of finite square wells, as illustrated in Fig. [2J On the other hand, the 
"heteroj unctions" between the molecular frame and the left and right vacua can be represented 
by infinitely high abrupt potential barriers with good approximation |101 \TT\ [T2] (not shown 
in Fig. [2]). Since the carriers are confined in the three directions, this heterostructure is zero- 
dimensional and should be called 'quantum dot' instead of 'quantum well'. 

Now, in each region the electron and hole quantum states are governed simply by corresponding 
one-dimensional effective-mass Schrodinger equations [31 [T21 [131 [2] 



n 1 



d 2 



2m* (x) dx 2 



+ Vi(x) 



ipi(x) = Si^i(x) 



(1) 



CD 
LU 



-8.7 



-9.0 



-9.3 



Figure 2. Top: molecular structure of the quantum dot PPP8-PDA7-PPP8. Bottom: Electron potential 
profile and energy levels calculated by means of the Extended Hiickel method. The barriers associated with the 
"heterojunctions" between the molecular frame and the left and right vacua are not shown. 



where x is the coordinate along the longitudinal direction of the molecule, and Vi(x), m*{x), 
and ipi(x) are the potential energy, (position-dependent) parabolic effective mass, eigenenergy 
and EMO of the carrier (i = electron or hole), respectively. For electrons (holes), the height of 
the barriers is given by the difference, or 'offset', between the energies of the LUCO's (HOCO's) 
of the parent polymers. The EMO approximation comes in handy because the bound electron 
and hole states can be estimated by the familiar particle-in-a-box expression |13t HH] 



7r 2 h 2 nf 

Ei~Ei + -f , (2) 
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where e« is measured from the edge of its respective band, Ei, m* is the carrier parabolic effective 
mass of PDA and l w is the width of the potential well. 



That this heterostructure is of type I means that the potential wells for electrons and holes 
are aligned ('straddling'). This implies that the superlattice can exhibit mini-bands both for 
electrons and holes. This report considers electron mini-bands only, the hole mini-bands being 
analogous. Moreover, attention will be focused on the mini-bands arising from two parent elec- 
tron levels bound to the intramolecular well. With the help of Eq. ([2]), it was predicted that the 
quantum dot PPPs-PDAy-PPPg contains two electron bound states, which were subsequently 
calculated more accurately employing the quantum-chemical methodology indicated in Section 
[3l The chemical structure, electron potential profile and energy levels of this quantum dot are 
shown in Fig. [2j The length of the molecule is 10.2 nm. 

Before proceeding to the superlattice, it is instructive to examine the evolution of the energy 
levels as the number of coupled quantum dots increases. Fig. [3] illustrates such an evolution for 
periodic arrays of quantum dots with molecular structures PPPs(-PDA7-PPPg) n , determined 



employing the methodology indicated in Section[3l As expected, it is observed that a manifold of 
n levels emerges from each parent level. In addition it is seen that for n = 15 the two manifolds 
already appear as quasi-continuous mini-bands. Naturally, in the hypothetical limit n — > oo the 
manifolds will become continuous mini-bands [3"1 I13| 
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Figure 3. Evolution of the energy levels in the periodic arrays of increasing lengths PPP8(-PDA7-PPP8) n 
calculated by means of the Extended Hiickel method. For reference, the levels of each array are drawn inside the 
potential profile of the parent quantum dot. 



Within the framework of the EMO theory, step (3), the assembly of the superlattice from the 
quantum dot just designed, is straightforward. The effective potential profile for electrons in 
the superlattice takes the form of a Kronig- Penney model [2], i.e. a periodic array of square 
wells, as displayed in Fig. [H It is evident that two parameters, the length of the well (l w ) and 
the length of the barrier (lb), are available for mini-bandstructure tailoring purposes. l w largely 
determines the energies of the parent discrete levels, allowing to perform 'level engineering'. 
Consequently, this parameter will also control the locations of the mini-bands and partly the 
(mini-)gap between them, permitting to carry out 'bandgap engineering'. On the other hand, 
lb determines the size of the overlap between adjacent parent EMO's. Hence, it will control the 
mini-bandwidths, and, consequently, the mini-effective masses, allowing to perform 'effective- 
mass engineering'. 
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Figure 4. Potential profile for electrons in the superlattice (PDAt-PPPs)™- The energy is measured from the 
LUCO of PDA. 

The Kronig-Penney model already helps to explain qualitatively the behavior of the energy lev- 
els displayed in Fig. [3l Specifically, the first manifold is narrower (denser) than the second 
because the overlap between adjacent EMO's of the ground levels of the individual quantum 
dots is smaller than for the excited ones. This is so because the height of the barrier is effectively 
smaller for the excited levels, making the tunneling probability between them larger. 

The methodology for the determination of the mini-bandgaps and electron and hole mini-effective 
masses as functions of l w and lb, entailed in steps (4) and (5), will be explained in the next section. 

3. Computational Methodology 

The bandstructures of the parent polymers PPP and PDA were calculated employing the 
semiempirical methodology of the previous contribution to this issue [8] , which entails a com- 
bination of the Austin Model 1 [15j . for geometry optimizations, with the Extended Hiickel 
method [16], for single-point electronic-structure calculations. Such methodology yields accu- 
rate HOCO and LUCO energies, but does not guarantee accurate bandwidths or, consequently, 
effective masses. Therefore, the PPP and PDA electron and hole parabolic effective masses were 
determined by means of an accurate real-space technique, based on interpolations of HOMO and 
LUMO data for oligomers of increasing sizes, which is explained elsewhere [10L fTTl [T2] . The elec- 
tronic structures of the PPP/PDA/PPP quantum dot and the finite superlattice PPPs(-PDA7- 
PPPg)^ were calculated using a methodology analogous to the one employed for copolymers in 
ref. [8], which yields somewhat less accurate HOMO and LUMO values than the HOCO and 
LUCO ones produced by the above-mentioned method employed for the parent polymers. 

The mini-bandstructure of the superlattice (PDA x -PPP y ) n can be determined employing pe- 
riodic boundary conditions and any quantum-chemical electronic structure method (see the 
previous paper of this issue [S] and references therein). However, for construction of calibration 
curves as functions of x and y, this would be very tedious and time-consuming. Furthermore, 
these methods produce a lot of atomistic information that is not required for the present pur- 
poses. On the other hand, the EMO method, employing the Kronig-Penney model, permits the 
calculation of the mini-bandstructure with reasonable accuracy, as long as the blocks are not 
too short [12], and very little computational effort [10 } lll[ [T2], as is now explained. 

In a finite superlattice, the carriers tend to be more localized towards the center of the structure, 
although states localized at the chain ends (Tamm states) may appear [14] . On the other hand, 
in an infinite superlattice the carriers are always equally likely to be found in any of the wells. 



Consequently, the envelope orbitals are periodic and satisfy Bloch's theorem [3j [T3] 

ip(x,k) = u(x,k)exjp(ikx). (3) 

where u(x, k) is a background function that possesses the same periodicity as the potential and k 
is the wavenumber along the longitudinal direction x. The first mini-Brillouin zone corresponds 
to —tt/L < k < tt/L, where L = l w + Z& is the period of the potential. With the exception 
of the possible appearance of Tamm states, an infinite superlattice is an adequate model for a 
long, finite superlattice. Thus, for simplicity, that model will be considered from this point on, 
and its envelope functions will be called "envelope crystalline orbitals" (ECO's) to differentiate 
them from the localized EMO's. 



Within the well and barrier regions the ECO's can be written as [13] 

ifi w = Aexp(ik w x) + B exp(— ik w x) (4) 
ipb = C exp(ikbx) + D exp(— ikbx) (5) 

respectively, where k w = ^2m^ ] e/h 2 and kb = \j2m* h {e — V)/H 2 , with m* (r = w,b) the effec- 
tive mass of the charge carrier inside the region r. 

To guarantee good behavior of the ECO's, these expressions must be matched across the two 
PPP-PDA heterojunctions present in the unit cell, enforcing the BenDaniel-Duke boundary 
conditions O [THJ HU [T7] that both ip(x) and be continuous at these points. Moreover, 

the periodicity condition [TT] 

ip(x + L) = ip(x) exp(ikL). (6) 

implied by Bloch's theorem (Eq. ([3])) must also be taken into account. (Notice that the Bloch 
wave number k gives the phase change of the ECO's from one unit cell to the next, whereas 
the wave number k r (r = w,b) governs the variations of these functions within each region). 
Employing these conditions, four coupled homogeneous equations are obtained, which, after 
some algebra, yield the transcendental equation [T4] 

— ^— — I = cos(kL) (7) 

2m w mbK w K J 

where k = ^2m h *(V — e)/h 2 . The values of the effective masses that appear in this equation are 
taken from ref. [TO] and are displayed in Table [TJ The superlattice dispersion curves e(k) were 
obtained from graphical solutions of Eq. ([7|). 



4. Results and Discussion 

Fig. [5] displays the two lowest-energy solutions of Eq. ([7]), in the extended mini-zone scheme, for 
the prototypical superlattice (PPPs-PDA7) n _ i>00 . It is worth noticing at the outset that the top 
and bottom of the first and second mini-bands, respectively, occur at the edge of the mini-zone. 
This implies that this kind of structure can be used for the design of inter-mini-band lasers, 
which are high-power sources of mid-infrared radiation |18j . 

Analogously to bulk crystals, a gap between two mini-bands is defined as the difference between 
the top of the lower-lying one and the bottom of the higher-lying one. By the same token, 
a mini-bandwidth is defined as A = \e(k = 0) — e(k = ir/L)\. The calculated mini-bandgap 
and the bandwidths of the first and second mini-bands of this superlattice are presented in 
Table [TJ The corresponding results of the electronic-structure calculation for the small super- 
lattice PPPg(-PDA7-PPPg)i5 shown in Fig. [3] are also presented in Table [TJ These two sets 




Figure 5. Mini-bandstructure of the infinite superlattice (PPP8-PDA7) n _>oo, calculated employing the ECO 
approximation and the Kronig-Penney model. The energy is measured from the LUCO of PDA. 



of results should not be compared directly, since, as mentioned in Section [3l the semiempirical 
electronic-structure methodology employed for this heterostructure produces somewhat inaccu- 
rate HOMO-LUMO gaps and unreliable bandwidths (and, consequently, effective masses). That 
is probably why the mini-bandgap of PPP8(-PDA7-PPPg)i5 turns out to be smaller than the 
one of PPP8(-PDA7-PPPg) n _ s . 00 by about 30%, and the first and second mini-bandwidths of the 
former are off by factors of three and four, respectively, with respect to the ones of the latter. 
However, the semi-quantitative agreement between these results is indicative of the plausibility 
of the ECO method. For comparison, the bandgaps and conduction-bandwidths of the parent 
polymers are also reported in Table [TJ It is observed that the superlattice mini-bandgaps fall in 
the mid-IR, whereas the parent polymer bandgaps fall in the UV-Vis. 

Within the parabolic approximation, the effective mass of a particle with a dispersion relation 
e(k) can be obtained from ^ = -^p^ [31 fT3"J I14j . The mini-effective masses of the electron 
around k = 0, relevant for direct transitions, for the first and second mini-bands, extracted by 
means of this equation from the dispersions of Fig. [5J are reported in Table [TJ It is interesting 
to compare these values with the electron effective masses of PDA and PPP (see Table UJ) |lUj . 
It is seen that the effective mass of the first mini-band is much larger than the "bulk" values, 
whereas the absolute value of the effective mass of the second mini-band is less than the bulk 
values. Since a carrier's mobility is inversely proportional to its effective mass [3], this implies 
that electrons in the second mini-band will have a greater response, and in the opposite sense, 
to an electric field than in the bulk, whereas the contrary occurs for the first mini-band. From 
now on, for the effective mass of the second mini-band only the absolute value will be reported. 



Table 1. Bandgaps (E g ), conduction-band widths (A2) and electron effective masses (m^) of 
the parent polymers; and mini-bandgaps (E g ), mini-bandwidths (Ai, A2) and mini-effective 
masses {m\,m*£) of the prototypical finite and infinite superlattices. mo is the electron mass. 
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(PPPg-PDAr)^ 


796 


4.11 


117 


0.818 


-0.015 



Let us now examine the effect of increasing the number of monomers in the barrier (PPP) 
regions over the mini-bandstructure for the superlattices (PPP :r -PDA7) ra _ s . 0O , with x = 5, .., 10, 
according to the ECO method employing the Kronig- Penney model. Fig. [6] shows the behavior 
of the extrema of each mini-band, where the lines joining the data, as in all the forthcoming 
figures, are obtained by interpolation. 
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Figure 6. Allowed energies (region between the same kind of line) of the infinite superlattices (PPP^- 
PDAt^-j-oo, where x was varied from 5 to 10. The energy is measured from the LUCO of PDA. The oligomer 
PPP8 corresponds to a barrier width of ~ 3.5 nm. 



As expected, it is seen that these get closer together as the barrier thickness increases. For a very 
thick barrier all the states within the mini-band will finally merge towards the corresponding 
level of the single well. This is a crossover from a superlattice to an array of non-interacting 
identical quantum dots. The left part of Fig. [7] displays the behavior of the mini-bandwiths, 
where it is observed that the convergence toward zero is monotonic. 
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Figure 7. Infinite superlattices (PPP a; -PDA7)„_» o. Left: mini-bandwidths. Right: mini-effective masses in 
units of the PDA electron effective mass (dotted line). 



The right part of Fig. [7] shows the behavior of the mini-effective masses, where it can be ap- 
preciated that, as compared to the one of the first mini-band, the effective mass of the second 
mini-band is largely insensitive to the barrier thickness, remaining below the bulk PDA value. 
This behavior is due to the proximity of the second mini-band to the edge of the barrier, which 
endows these states with a much stronger PPP character than the states of the first mini-band, 
making them much less sensitive to variations in the thickness of the barrier. 

Let us finally examine the effect of increasing the number of monomers in the well (PDA) regions 
over the mini-bandstructure for the superlattices (PPPg-PDAy),^,^, with y = 5, .., 10. The left 
part of Fig. [8] displays the results for the mini-bandwidths, where it is observed that these 
decrease as the width of the well grows. This can be explained by the fact that as the width of a 
single quantum dot increases the energies of its confined states decrease and the corresponding 
EMO's become more compact, thereby decreasing the overlap between adjacent levels in the 
superlattice. 



The right part of Fig. [8] presents the results for the mini-effective masses. Again, it is seen 
that the effective mass of the second mini-band is largely insensitive to the barrier thickness. 
However, now it begins to grow at about 6.3 nm and surpasses the PDA bulk value at about 
7.3 nm. This can be understood by realizing that as the well gets wider the states of the second 
mini-band steadily move further down from the edge of the barrier, steadily losing their PPP 
character and gaining more PDA character, which makes them sensitive to the well width. This 
also explains the appearance of a shoulder in the mini-bandwidth curve around 6.3 nm, which 
must be due to the crossover from dominating PPP character to dominating PDA character. 

The ECO approximation breaks down in the vicinity of abrupt heteroj unctions |12j . Therefore, 
the more similar the chemical structures of the parent polymers are and the longer the blocks 
of the superlattice become, so that the effect of the heteroj unctions is effectively diminished, 
the better this approximation works [12J. To assess the accuracy of the ECO method more 
quantitatively, a comparison against higher-level bandstructure calculations is necessary. 




Figure 8. Infinite superlattice (PPPs-PDA y )„. Left: mini-bandwidths. Right: mini-effective masses in units 
of the PDA electron effective mass (dotted line). 



5. Summary, Conclusions and Perspectives 

With the aid of the envelope molecular-orbital theory, molecular organic quantum-dot super- 
lattices with the general structure (PDA x -PPPj / ) n were designed specifically to display two 
mini-bands with energies in between the conduction-band edges of the parent PPP and PDA 
polymers. Calibration curves of the gap between these two electron mini-bands and their effec- 
tive masses as functions of the dot and barrier lengths were constructed, which allow tailoring of 
these heterostructures for specific electronic or optoelectronic applications, like infrared detectors 
and lasers. The mini-bandstructures were determined employing an envelope crystalline-orbital 
approximation with a Kronig-Penny model, which is seen to be a conceptually appealing and 
computationally very economical methodology. 

As expected, it was found that the lower-lying mini-band exhibits a smaller bandwidth than the 
higher-lying one, and that these bandwidths decrease as the lengths of the dot and the barrier 
increase. Perhaps not so expectedly, it was also observed that the effective mass of the higher- 
lying mini-band is largely insensitive to increments in the dot and barrier lengths, remaining 
below the bulk PDA effective mass in most cases, whereas the effective mass of the lower-lying 
mini-band increases rapidly with such increments. 

In future studies a quantitative comparison between the crystalline-orbital approximation and 
higher-level atomistic bandstructure methods will be performed. In addition, the full mini- 
bandstructure, i.e. including the hole minibands, will be considered, so that interband transitions 
and mini-excitonic effects in these molecular organic heterostructures can be determined. 
Furthermore, the effects of negative and positive doping of the superlattices, i.e. reduction and 
oxidation of the molecules, respectively, on their electronic, optical and conduction properties 
will be investigated. 
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